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Abstract. We present a general and systematic theory of non-equilibrium dynamics of multi-component 
fluid membranes, in general, and membranes containing transmembrane proteins, in particular. Developed 
based on a minimal number of principles of statistical physics and designed to be a meso/macroscopic- 
scale effective description, the theory is formulated in terms of a set of equations of hydrodynamics and 
linear constitutive relations. As a particular emphasis of the theory, the equations and the constitutive rela- 
tions address both the thermodynamic and the hydrodynamic consequences of the unconventional material 
characteristics of lipid-protein membranes and contain proposals as well as predictions which have not yet 
been made in already existed work on membrane hydrodynamics and which may have experimental rele- 
vance. The framework structure of the theory makes possible its applications to a range of non-equilibrium 
phenomena in a range of membrane systems, as discussions in the paper of a few limit cases demonstrate. 
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1 Introduction 

Lipid-protein fluid membranes are the most essential struc- 
tural element in biological cells, defining boundaries of 
the cells and intracellular organelles. They are also one of 
the important functional elements, actively participating 
in many cellular processes. Each biological membrane is 
composed of a core bilayer of amphiphilic lipids, in which 
transmembrane proteins are embedded and with which 
peripheral proteins are associated. In its functional state, 
the membrane is fluid, allowing the constant movement 
and organization of the constituent molecules within the 
structure, and its two-dimensional geometry deforms eas- 
ily in connection with its function, requiring energies com- 
parable to thermal energy only. Moreover, the membrane 
constantly exchanges material and energy with its envi- 
ronment; active transport of small solute molecules across 
the membrane takes place constantly, carried out by mem- 
brane proteins that require external energy sources, such 
as chemical energy provided by ATP, electrochemical en- 
ergy stored in cross- membrane proton gradients, or light. 
From the point of view of statistical physics, it is obvious 
that such a functioning membrane should be treated as 
a non-equilibrium system and that the dynamics of the 
membrane is intimately coupled to the dynamics of the 
bulk fluids within which it is embedded. 
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During the last three decades, investigations, charac- 
terizations and understanding of the physical properties 
and behaviour of such membranes and simpler model mem- 
branes have fueled both the development of statistical 
physics of soft condensed matter in general and the de- 
velopment of membrane biophysics in particular, because 
these systems have posed many issues challenging the tra- 
ditional framework of statistical physics and also because 
it has become more and more appreciated that under- 
standing of the physics of the membranes can shed light 
on their biological functions 

For obvious reasons, most of the development of mem- 
brane statistical physics has focussed only on the equi- 
librium, static aspect of the membrane systems. To the 
purpose of a more complete physical description of biolog- 
ical membranes, however, their non-equilibrium behaviour 
must be investigated and described. The work presented 
in this paper is our attempt at taking a step in that di- 
rection. 

Some recent experimental as well as theoretical stud- 
ies of simpler model systems of biological membranes have 
particularly motivated our work ^ . In the model systems, 
a single type of transmembrane protein, which can be ex- 
ternally driven to actively transport small ions across the 
membrane, was reconstituted into a core lipid bilayer at 
various concentrations. By turning on or off the exter- 
nal driving force, the lipid-protein composite membrane 
was then set in either a non-equilibrium or an equilibrium 
state, and the non-equilibrium dynamics of the membrane 
conformation was then investigated experimentally. The- 
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ories, in the form of equations of dynamics, were also de- 
veloped to describe or interpret the experimental observa- 
tions. The formulation of the theories appears, however, 
entirely intuition based, which makes its validity rather 
opaque and its generalization to other lipid-protein com- 
posite membranes difficult. We have, therefore, made it 
our goal to develop a more general and systematic theo- 
retical description of lipid-protein composite membranes 
in non-equilibrium states, believing that such a theoretical 
framework will become useful as more and more experi- 
mental studies of such systems will be done with their 
results to be interpreted theoretically. 

In fact, a conceptual basis for developing a general de- 
scription already exists, and it consists of a minimal num- 
ber of principles of statistical physics: the basic laws of 
thermodynamics, formulated as a balance equation of en- 
tropy, the assumption of local thermodynamic equihbrium 
in a non-equilibrium system and the dynamic formulation 
of laws of conservation of mass, momentum and energy, 
and if relevant, angular momentum.^ However, the ex- 
plicit expressions and implementation of these principles 
in the case of lipid-protein composite membranes, which in 
the end result in a set of clearly formulated equations of 
(hydro)dynamics, must be developed systematically and 
unambiguously to reflect the physical characteristics per- 
taining to the membranes. 

At the level of continuum descriptions, a membrane is 
modelled effectively as an interface separating two bulk 
fluids. Its physical characteristics, however, distinguish it 
from a simple interface separating two coexisting fluids. 
Firstly, its molecular composition and material structure 
differ significantly from those of the two bulk fluids with 
which it is in contact, and the motion and organization 
of its constituent molecules within its interface structure 
provide a set of interesting phenomena. Thus, the material 
content, and associated with it, the energy content as well 
as momentum content of the membrane should in general 
be taken into account. Secondly, its energy (or, thermo- 
dynamics) depends not only on its surface area, but also 
on its local surface geometry, or curvatures 3 4 . Thirdly, 
its properties with regard to processes of transport of ma- 
terial and energy are in general markedly different from 
those of the bulk fluids. A good example is the much 
slower diffusion across the membrane than in the bulk flu- 
ids of polar molecules. Thus, the dissipations associated 
with such transport processes should not be neglected a 
priori, as is canonically done when non-equilibrium dy- 
namics involving a conventional interface is considered. In 
fact, there has been experimental evidence suggesting that 
a membrane should not be considered as dissipationless in 
its non-equilibrium state j5l()l7j . Finally, the typical con- 
stituent molecules of a lipid-protein composite membrane 
- the smaller amphiphilic lipid molecules and much larger 
transmembrane protein molecules - differ significantly in 
their molecular chemistry and structures. Consequently, 
their interactions with the contacting bulk fluids may also 

^ These conservation laws are sufHcient if we limit ourselves 
to the non-equilibrium dynamics which does not involve non- 
conserved physical variables. 



be expected to differ. Our work, in essence, consists in 
recognizing these unconventional characteristics, clarify- 
ing the basic conceptual issues that inevitably rise from 
considering them, and finally developing an unambiguous 
description of these characteristics in the context of non- 
equilibrium dynamics of the system. 

To be sure, studies concerning fundamental descrip- 
tions of non-equilibrium dynamics of fiuid surfaces and 
interfaces began already three decades ago |8I9| . For ex- 
ample Bedeaux et al. developed for conventional inter- 
faces a formulation of non-equilibrium dynamics by taking 
into account the energy and momentum content of the in- 
terface, from which a connection to equilibrium interface 
thermodynamics can be made [H] ■ More recently, a theory 
for the non-equilibrium dynamics of a two-component sur- 
factant interface separating air and water was presented 
in the form of a set of dynamic equations, with an em- 
phasis on nonlinear phenomena where phase separation 
and surface deformations are coupled In this the- 
ory, the mass content and the composition of the inter- 
face were explicitly taken into account, and the depen- 
dence of the interface thermodynamics on the curvatures 
was also introduced. The derivation of the dynamic equa- 
tions was, however, more intuitive than systematic, and 
certain assumptions were made, which were neither nec- 
essary nor justified. Moreover, the transverse transport 
processes were not discussed. A generalization of this the- 
ory to lipid-protein composite membranes did not appear 
straightforward. 

It should be helpful to the reader that we at the outset 
describe more specifically the system and its non-equilibri- 
um state that we have in mind when developing our the- 
ory. The essential picture is sketched out in Fig. 1(a). The 
system we consider consists of a multi-component mem- 
brane of lipid-protein composite, which has a nanometer 
thickness and which assumes a quasi two-dimensional geo- 
metrical shape, and of two aqueous fiuids, which the mem- 
brane is in contact with. The two fluids can be identical 
in their chemical compositions and equilibrium thermo- 
dynamics, or distinctly different, as in the case of two 
fluids under conditions of phase separation. In order to 
keep the theory general, we allow the membrane geome- 
try to be arbitrary, and we do not specify quantitatively 
the number of molecular species constituting the mem- 
brane, nor their exact chemical structures, although the 
cases where the protein components are transmembrane 
proteins are of particular interest to us. We also assume 
that the aqueous fluids may contain more than one types 
of solute molecules. 

We then consider a general non-equilibrium state, where 
chemical, thermal and mechanical gradients exist both in 
the directions tangent to the membrane surface and in the 
direction transverse to the membrane, and where chemi- 
cal reactions may take place on the membrane. But in the 
scope of this paper, we will limit ourselves to situations 
where the length scales over which the assumption of local 
thermodynamic equilibrium is valid are larger, though not 
by orders of magnitude, than the thickness of the mem- 
brane. The corresponding time scales may be expected to 
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Fig. 1. A schematic sketch of the membrane-fluid system un- 
der our consideration. Fig. 1(a) depicts the interfacial region 
in the real system, showing a membrane composed of bilayer- 
forming lipids and transmembrane proteins in contact with two 
fluids containing small solutes. Fig. 1(b) illustrates the repre- 
sentation of the interfacial region in the corresponding Gibbs 
model system. 

approach the mesoscopic or above. Given the typical time 
resolutions of around milliseconds of most experimental 
techniques used to study membranes, our theory will have 
a reasonably large range of applicability in terms of length 
and time scales. It follows then that the membrane is as- 
sumed to be in thermodynamic equilibrium with the fluids 
which it is in immediate contact with.^ 

We would also like to briefly state the philosophy which 
we have followed when developing the theory. Although 
the membrane-fluid systems we consider are on micro- 
scopic length scales highly inhomogeneous in molecular 
distributions and transport properties, our theory is not 
concerned with accurately describing the details of the 
microscopic-scale inhomogeneity, but aims at describing 
the inhomogeneity only through those effects that may 
be observed at mesoscopic or macroscopic scales. Follow- 
ing this philosophy, we have employed an idea originated 
from Gibbs ^J: the idea of replacing a real, inhomoge- 
neous membrane-fluid system with a model system of two 
homogeneous bulk fluids separated by an infinitely thin di- 
viding surface and relating to the dividing surface all the 
excess thermodynamic and hydrodynamic contributions 
that are required by an equivalence between the thermo- 
dynamic and hydrodynamic behaviour of the two systems. 
This idea is also illustrated graphically by Fig. 1. For later 
convenience, we will dub this model system "the Gibbs 
system." 

The paper is organized as follows. In Section |2 the 
necessary surface differential geometry is briefly reviewed. 
In Section |31 the concept of the surface thermodynamics 
of the membrane is unambiguously defined based on the 
idea of the Gibbs model system, and a general expression 
of the surface thermodynamics is derived and discussed. 
In addition, a general identity relating different thermo- 
dynamic variables associated with the membrane surface 
is derived, which takes into account the salient mechani- 
cal characteristics of membranes. To our knowledge, this 

^ Formulation of thermodynamics in situations where this 
assumption is not valid was first discussed in Ref. [S| in the 
context of "dynamic surface tension." 



is the first derivation of such an identity. An assumption 
made in this Section is that intrinsic orientational degrees 
of freedom are not relevant. In Section^ a description of 
the hydrodynamics of the whole membrane-fluid system 
is formulated by deflning relevant "bulk" and "surface" 
hydrodynamic variables and by establishing equations of 
dynamics for all of the hydrodynamic variables based on 
the fundamental conservation laws. The hydrodynamic de- 
scription is limited to cases, where systems under consid- 
eration have no intrinsic angular momenta. In Sectional 
the results from Section |31 and Section 01 combined with 
the assumption of local thermodynamic equilibrium are 
used to derive the entropy production, based on which 
conjugate pairs of general thermodynamic/hydrodynamic 
force and dissipative current arc identified. Constitutive 
relations for the linear non-equilibrium dynamics of the 
system are derived in Section |H| In particular, some of 
these constitutive relations, which under appropriate con- 
ditions reduce to those that appear familiar, make predic- 
tions about new mechanisms governing various dissipative 
processes in the membrane. The complete set of consti- 
tutive relations, together with the equations of dynamics, 
thus provide a closed formulation of the hydrodynamics of 
the whole membrane-fluid system. In Section [3 a number 
of limit cases of the general theory, which have practical 
relevance, are discussed, in order that the theory, having 
been presented in a general and formulistic way, also be 
seen from a practical point of view. Finally, in the con- 
cluding section. Section |H1 the theory is discussed within 
the context of its applications and its connections to ex- 
perimental measurements. In order that a technical point 
be made clear, which we expect will often be encountered 
in applications of the theory, a short appendix is also at- 
tached. 



2 DifFerential geometry of surfaces 

In the continuum theory which we will develop in this 
paper, a membrane is structurally described as a two- 
dimensional surface. In this section we briefly review, main- 
ly to establish the notation, the mathematical language of 
two-dimensional differential geometry, which will be used 
to describe membrane geometry. A more comprehensive 
introduction can be found in, for instance, Refs. jl2ll[{| . 

The dynamic shape of the surface is represented by 
a space- vector function R — il(^^, t). The variables 
and are internal coordinates corresponding to a parame- 
trization of the surface and t represents time. At each 
point on the membrane a basis for three-dimensional vec- 
tors can be established. Two of them are the tangential 
vectors 

t^^d^R^^, (1) 

where a = 1,2, and the third is a unit vector normal to 
the surface, 

^ diR X d2R , . 

" ~ \diR X 92^1 ■ ^ ' 
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The local geometry of the surface is characterized by 
two surface tensors, the metric tensor and the curvature 
tensor. The local metric tensor is defined by 

9afi = daR ■ dpR . (3) 

It has an inverse, g"^ , which satisfies by definition 

9"^9f^, = , (4) 

where 5" is the Kronecker delta and where the repeated 
Greek superscript-subscript indices imply summation fol- 
lowing the Einstein summation convention. The metric 
tensor and its inverse are used to raise and lower Greek 
indices as in the following example: 

= g'^^tp , = go^pt^ . (5) 

The curvature tensor Kap is defined via the second deriva- 
tives of the surface shape function, 

Kap = n ■ dadpR . (6) 

From Kaf3 the scalar mean curvature H and Gaussian 
curvature K can be obtained: 

H = \g^^K^p , (7) 

K = det g°''^Kp^ . (8) 

Two other tensors will also be introduced here for later 
convenience, 

= ea/3V^ , e"^ = e'^^l^ , (9) 

where = e"^ with en — €22 — and ei2 = —£21 = 1 
are relative tensors, and g = det g^p is the determinant of 
the metric tensor. 

Expressions of covariant/contravariant differentiations 
of vector and tensor functions defined on the surfaces are 
facilitated by the use of the Christoffel symbols, r^p- One 
instance, which will become particularly useful later, is 
the covariant differentiation of a surface vector function, 
w = w°'ta, given by 

Do^w^ = da^wP + wT^^ . (10) 

The Christoffel symbols can also be defined as certain 
combinations of the derivatives of the metric tensor, 
namely, 

^7 _ 1 7^ ( dgsa dgps dg^p \ 

It follows that the covariant divergence of can be writ- 
ten as 

Dcw^ = . (12) 

Finally, the area of a local differential element of the 
surface is given by 

dA = ^dede , (13) 

an expression which will be repeatedly used in surface 
integrals. 



3 Membrane thermodynamics 

An important component of a general description of non- 
equilibrium dynamics of the membrane-fiuid system sketch- 
ed out in Introduction is an appropriate formulation of 
the equilibrium thermodynamics of the system. In this 
section we will present such a formulation. Two points 
which pertain to the system and which have been dealt 
with in the formulation, may already be mentioned here. 
Firstly, to describe the thermodynamic effects arising from 
the microscopic inhomogeneity inherent in the system, we 
have employed the idea of the Gibbs model system. A sub- 
tle issue of principle arises, however, when Gibbs' idea is 
applied to the membrane system. In the case of a con- 
ventional capillary interface, whose mechanical property 
is entirely described macroscopically by a surface tension 
experimentally measurable, the position of the Gibbs di- 
viding surface is uniquely determined by the thermody- 
namic - including mechanical - equivalence between the 
real and the Gibbs systems 14 15 . In the case of the mem- 
brane system, whose macroscopic mechanical properties 
include not only a "membrane tension," but also "mem- 
brane bending moments," the thermodynamic equivalence 
between the real and the Gibbs systems still leaves the po- 
sition of the Gibbs dividing surface free to be chosen in 
principle. With a view of making connection to experimen- 
tal, in particular mechanical, characterizations of mem- 
brane systems, where the geometrical profile of a mem- 
brane surface is resolved with optical resolutions we 
do not define quantitatively the position of the theoret- 
ical surface, but will work under the assumption that a 
position can be chosen to be consistent with the experi- 
mentally determined one. 

Secondly, related to the significance of membrane bend- 
ing moments in the effective description of a general mem- 
brane system | 3I4| , the thermodynamic free energies of the 
membrane systems under our consideration must depend 
on local geometrical properties such as the principal curva- 
tures. We recognize that, due to such dependence, the free 
energies no longer scale homogeneously with the size of the 
membrane and that, consequently, the Gibbs adsorption 
equation in its canonical form |14j no longer exists. How- 
ever, a certain statement, in the form of an identity, can 
still be made about how the different surface thermody- 
namic variables are related. The derivation of this identity 
will also be given in this section. 

3.1 The basic equation 

In order to develop a local description of the equilibrium 
thermodynamics of the membrane system, it is necessary 
that we consider a small cell, or a volume element, i7, of 
the whole membrane-fluid system. A sketch of the cell is 
given in Fig. 2. The spatial extensions of the cell are de- 
fined by a base area element S on the dividing surface, 
which spans in a chosen internal coordinate space a fixed 
rectangle with corners at (^^ - A£_^/2,£,'^ ± and 
+ Z\^V2,C^ ± ^C^/2), and by a height e+ above and 
a height t~ below the dividing surface in the direction of 
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Fig. 2. An illustration of the cell, or volume element, S. 
and refer to the side surfaces of the cell. 



the surface normal vector n. The local geometry of the 
sm'face is described by e+ and e~ are chosen 

such that the physical characteristics of the system at e+ 
and e~ reach those of the two homogeneous bulk fluids. 
The cell is assumed to be in thermal equilibrium with a 
uniform temperature, and it is also assumed to be in me- 
chanical equilibrium, although this does not necessarily 
imply a uniform pressure across the whole cell due to the 
fact that the membrane interface may be curved. More- 
over, in order to include mass motion in the formulation 
of thermodynamics, the cell is considered to be in uniform 
motion characterized by a velocity v. 

Regarding the state of chemical equilibrium in the cell, 
particular considerations are needed. It is well known that 
a typical membrane appears impermeable to polar mole- 
cules such as small ions on the time scales of seconds [T7| 
and that the reorientations of the constituent amphiphilic 
lipid and protein molecules from one side of the membrane 
surface to the other are not observed on similar time scales 
|18| . To take into account this fact, therefore, we introduce 
an operational definition of species, which is broader than 
the canonical definition based on the chemical nature of 
a molecular component: solution molecules of the same 
chemical structure, but found to be located on the two 
different sides of the dividing surface will be counted as 
two different species if their transport across the mem- 
brane is a slow process; lipid, or protein molecules having 
the same chemistry, but oriented in the two opposite di- 
rections across the membrane will also be considered as 
two different species. The equilibrium state which we will 
develop a thermodynamic description for corresponds to 
a state where each species is in its chemical equilibrium, 
but where there is no chemical equilibrium, in the quasi- 
static sense, between any two species of the same molecu- 
lar chemistry. To maintain the generality of the theory, we 
will also allow for in the system the presence of molecular 



components which do reach chemical equilibrium across 
the membrane in the same quasi-static sense. Each such 
component constitutes a single species by our definition. 

We denote the energy content in this cell by E^. Fol- 
lowing the formal structure of thermodynamics, we as- 
sume that E£ depends on only a few state variables: the 
molecular numbers of the different species, TV^ where 
"A" is an index labelling the species; the entropy S^; the 
momentum P^; and finally, the variables characterizing 
the shape and size of the cell, namely, the heights e^, 
which are assumed to be fixed, and derivatives of the shape 
field (such as H and K).^ The reason that only 

derivatives of R, not R itself, are allowed is that the ther- 
modynamics should be invariant under rigid translations, 
provided that the translational-symmetry-breaking effects 
such as the gravitational force are taken into account as 
external forces. Based on this assumption, we can there- 
fore express, for the chosen cell, the first and second laws 
of thermodynamics as a complete differential of with 
respect to all the relevant state variables: 

SE^^v 6Ps + J2 i^^^^A^E + T5Ss 

A 

-Fs-SR + d^{S%-6R) , (14) 

where we have already identified the partial derivatives 
with respect to Pjj, Nj^£, and S£ with the velocity v, 
the chemical potential for species A, and the temperature 
T of the cell. 

In the last two terms, is a regular quantity and 
carries the significance of a physical force, whereas 
contains differential operators and is related to surface 
stresses. From a mathematical point view, their presence 
is not difficult to understand, as they represent the vari- 
ation in the energy function resulted from any variation 
in the shape field, i2(^^,^^), and in turn, variations in its 
derivatives. Seen from a more physical point of view, the 
two terms must describe the mechanical work done on the 
cell when the shape of the dividing surface is changed. To 
illustrate how their functional forms arise from their me- 
chanical origins is not a trivial issue, and is discussed in a 
separate paper |20| . However, mechanical interpretations 
of and S'^ will be made a bit later in the paper, clari- 
fying the reason for expressing the work in those particular 
functional forms. 

To reformulate Eq. 114|l by the use of the Gibbs model 
system, we introduce bulk volume densities for extensive 
quantities of the cell: , , and . These can be ex- 
pressed as functions of intensive thermodynamic variables 
V, T and /i"^'s and are assumed to represent the known 
thermodynamic behaviour of the homogeneous bulk fiuids 
in the following sense: 

(5e± = • 5p± + T(5s± + ^ /i^^nj , (15) 

A 

In making this assumption, we exclude from our considera- 
tions cases where intrinsic orientational degrees of freedom and 
their independent rotational motion are relevant. An extension 
which does include those effects is given in Ref. |19I2()| . 
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and 

e±=vp±+rs± + ^^^n±-p± , (16) 

A 

where denotes the bulk function of velocity, tempera- 
ture and chemical potentials that corresponds to pressure. 

We can now define the excess area densities of exten- 
sive quantities, which must be associated with the dividing 
surface according to the construction of the Gibbs model 
system: 



e = 


1 






(17) 


P = 


1 






(18) 


s = 


1 


{s^-vp+- 




(19) 




1 






(20) 



where is the area of the surface element S, and 
represents the volume of the cell part which is above/below 
the dividing surface, respectively. It follows that 

p = pv, (21) 

where p = m^nA is the excess mass density. 
Using these excess quantities, together with 

1/| - As (^e± T {e^fH + \{e^fK^ , (22) 
As = AeAeVg , (23) 

where we have assumed that the radii of curvature are 
bigger than , we finally arrive at the following reformu- 
lation of Eq. ((TUl . 

6 (V^e) = v-S i^jp) + TS (V?s) + (VgriA) 

A 

- Vff/r. • SR + ^gD„ (S^ ■ 6R) . (24) 
Vector quantities f^^ and S" are defined by 

= ^ - D^Bl - D^BI , (25) 

5" = -^S% B" +n (cf + Cl'^^ dp , (26) 

where 

Bl= - e^p^t"' ± p±(e±)2 (^Hg°'P - \^'^^^ *'3 

-DpCfn, (27) 

Cf^ T{e^fp^\9-^ + \{e^fp^{2Hg'^P-K'^P) . 

(28) 

Note that the physics represented by Eq. 1)24(1 should be 
independent of any specific numerical values of e*. 



Based on Eq. ((24|l the mathematical area element S on 
the dividing surface may be viewed as an effective surface 
system which has its own thermodynamic properties and 
which interacts with its "surroundings." In this effective 
picture, the last two terms in Eq. - related to the 
mechanical-work terms in Eq. \1A\ - may be interpreted 
as the work done on the effective surface system by its 
surroundings: In Section 5, it will become clear that f^^ 
must balance the effective force exerted on S by sources 
external to the dividing surface under mechanical equilib- 
rium. Similarly, ^/gDa {S" ■ 6R) will be shown in Eq. 
to represent the work done on S by the rest of the divid- 
ing surface. But, this interpretation can already be made 
plausible here by noting that ^/gDa {S" ■ SR) can 

actually be rewritten as an integral over the boundary dS 
of the area element S 

[ (fS, y/^Da (S'" ■SR)= [ ds Va^S" ■ SR , (29) 

where s is the arc length and f^t" is a unit normal vector 
pointing away from the boundary dS. 

There is in fact a connection between /^^ and S"' . If 
5"o) is used to represent all the contributions in S" that do 
not involve differential operators, then it can be seen from 
Eq. (|24|l that the following relationship must be satisfied, 
as a consequence of the invariance of the thermodynamics 
under rigid translations: 

= i^^S^o) . (30) 

/j-s accounts for the total mechanical force exerted on the 
effective surface element S by the rest of the effective 
surface. In what follows, f^^ will be referred to as the 
"restoring force." It is clear that S"o) should be identified 
as the surface stress tensor as defined in |21I22| . 

The above interpretations make clear the reason for or- 
ganizing the geometry-dependent work explicitly into the 
two particular functional terms in Eq. 1(14(1 . We would like 
to note also that these interpretations can be obtained in 
a more physically intuitive and direct way by formulat- 
ing the mechanical work explicitly, once a model of the 
mechanical behaviour of the inhomogeneous cell is given. 
F£ and S"^, or /^.^ and 5", can be identified with the 
mechanical model quantities. But, we defer the discussion 
of that topic to another paper pUj . 

It can be seen easily that Eq. 1(24(1 does not define 
S°' uniquely. An addition to it of the following type, for 
instance, 

e''^ {dpV + Vdp) , (31) 

where V is an arbitrary vector, does not change the me- 
chanical work at all. This seemingly mathematical point 
implies in fact a non-trivial physical statement: the theo- 
retical characterization of the mechanical behaviour of a 
membrane system in terms of an effective surface stress 
tensor and surface bending-moment tensor is not unique, 
as opposed to the "belief implied in the canonical de- 
scription of membrane mechanics |23| . We will discuss and 
clarify this issue elsewhere (2()( . 

Eq. ((24(1 thus provides the basic equation of the mem- 
brane surface thermodynamics and will be used later in 
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our description of non-equilibrium dynamics of the mem- 
brane system. Once an explicit functional form of e and 
the values of the surface thermodynamic variables, p, s, 
tt-a's, and R are assumed to be known, Eq. (|^ . or its 
integrated form E = Jy^dA e, can be used to determine 
the other physically measurable thermodynamic variables 
of the surface as follows: 



V = 



T = 



frs 



1 SE 
1 SE 

1 SE 
1 SE 



R,{nA}-s 
p,R,{nA} 

p,s,R,{nB\B^A} 



(32) 



3.2 A useful identity derived from reparametrization 
invariance 

An important element of the canonical thermodynamics 
of a capillary interface between two coexisting fluids is 
the so-called Gibbs adsorption equation ^1], which relates 
together the different intensive thermodynamic variables 
defining the state of the interface. It results from the fact 
that the excess thermodynamic free energies associated 
with such an interface scale proportionally with the area 
of the interface at constant excess densities of the exten- 
sive quantities, or equivalently, that the mechanical char- 
acterization of the interface is given by a single intensive 
quantity, the surface tension. 

The thermodynamics of the type of membrane sys- 
tems under our considerations is, however, different. The 
canonical model of the membrane mechanics proposed by 
Helfrich and Evans provides an example. In the model, 
the part of the free energy associated with the membrane 
mechanics is given by dA (2kH^ + co) j where k and 
tTo SLie constants. It is clear that, if the linear size of the 
membrane surface is scaled by i2 — > XR, the term involv- 
ing CTo will increase with while the one including k (the 
bending term) will not change. Consequently, the Gibbs 
adsorption equation no longer exists; and the concept of 
surface tension alone is no longer sufficient to describe the 
mechanical behaviour of the membrane surface at meso- 
scopic or macroscopic scales. Instead, as Eq. (|24ll implies, 
the restoring force /^g provides an appropriate mechanical 
quantity. 

It is obvious from Eq. H24() that f^^ depends on the 
determining (surface) thermodynamic variables such as v, 
T, /i"^'s as well as the geometry of the dividing surface. It 
turns out that the tangential components of f^^ are inti- 
mately connected with the spatial inhomogeneities in v, T, 
and ^"*'s. This connection arises from the fact that both 
the thermodynamic and the hydrodynamic behaviour of 
the dividing surface is that of a two-dimensional "fluid sys- 
tem." In other words, they should remain invariant under 
any change of the internal coordinate system. 



To derive the explicit expression of the connection, 
we consider a situation where there exist spatial inhomo- 
geneities in the surface thermodynamic variables. The to- 
tal excess energy E associated with the dividing surface is 
then the surface integral of the local density e, i.e. a func- 
tional of p, s, UA^s and R. Under an arbitrary, infinitesi- 
mal change of internal coordinates, or "reparametrization" 
of the dividing surface. 



(33) 



the functional form of the local density of excess entropy, 
as an example of a physical quantity, must change from 
s(^^,^^) in the old coordinate system to another form 
s'iC ) in new such that 

s{e,e) - = s'{e,e) + scd^s , (34) 

as required by the reparametrization invariance. This is 
equivalent to making the following variation in the func- 
tional form of the entropy density expressed in the old 
coordinate system: 

ss EE s'{e,e) - sie^e) = -^rs^s . (35) 

The variations in the functional forms of p, R and ua can 
be expressed similarly. 

Under fixed boundary conditions, the above variations 
lead to a variation in the integrated energy E 

SE^ [ d'^ (v5{^p)+TS{^s) 
Jm \ 

- V^/,, -SR + Y, i^uA) ) , (36) 

where Eq. H24|) has been used. Obviously, this variation 
must be zero. 

Inserting into Eq. I|36|l Eq. (|35|l . its analogs for p, n^'s 
and ii, as well as i5(y^) = y^t" • daSR, and performing 
a partial integration yields 



SE=^ / dA [p-dav + sdaT 
Jm \ 

+ frs-tn+J2 '^^^"M^) Sr^O. (37) 

A ^ 

Since 5f " is arbitrary, it can finally be concluded that 
f,,-to,+p- dcv + sdcrT + ^ nA^,r^l^ = , (38) 

A 

must always hold. This identity will become a useful one 
in the formulation of non-equilibrium dynamics. 

Two remarks are worth making here. First, although 
the physical reason underlying the above identity appears 
conceptually obvious, we are not aware of any earlier work 
where the result has been systematically derived. Sec- 
ondly, although we have made the derivation with mem- 
brane systems in mind, the result also applies to a conven- 
tional capillary interface with a spatially varying surface 
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tension a, in which case fj.^ -t^ = da<J. Eq. I|38|l thus coin- 
cides with the expression of the Gibbs adsorption equation 
when it is appUed to cases where spatial inhomogeneities 
in the interface are present |14| . 



4 Dynamic Formulation of Conservation Laws 

Having formulated a description of local thermodynamic 
equilibrium of the membrane-fluid system in terms of the 
surface excess quantities, we now proceed to consider the 
general non-equilibrium state of the system as defined in 
Introduction and develop a theory which describes the dy- 
namics of the non-equilibrium state. Following the philos- 
ophy of developing an effective theory by the use of the 
Gibbs model system, where two bulk fluids meet at an 
infinitely thin dividing surface, we assume that the hy- 
drodynamic description of the two bulk fluids, in terms 
of their local thermodynamics and transport properties, 
is entirely known, and thus we relate, not only the ther- 
modynamic, but also the hydrodynamic, effects associated 
with the microscopic inhomogeneity in the real system to 
the dividing surface and to its modes of interaction with 
the two bulk fluids. 

To make the presentation easy to follow, we first define 
a few notations pertaining to the description of the space. 
Similar to previous notation, is used to repre- 

sent the dynamic shape of the dividing surface, and the 
space is then divided into two regions separated by the di- 
viding surface: -region refers to the bulk-fluid region 
which the normal of the surface, n(^^,^^,t) points into 
and "—"-region the other. A scalar function f{r, t) is intro- 
duced such that it is zero on the dividing surface and posi- 
tive/negative on the +/— side, respectively; two Heaviside 
step functions are then defined as 9'^{r,t) = 9{±f{r,t)). 
A few identities follow immediately, 

dtO^ = T / dAn-dtR6{r~R) , 
Jm 

V9^ = ± I dAnS{r-R) , 
Jm 



to, ■VS{r-R) = -dc5 {r-R) , 



(39) 



which will be used below. 

The starting point of the hydrodynamic description 
is a formulation of the basic laws of conservation of the 
molecular number of each species, momentum and energy 
for the model system in the context of non-equilibrium 
dynamics. Specifically, it is assumed that the following 
equation of dynamics holds. 



5 / N ^ 



Jx{r,t) + 



(40) 



where a;(r, t) represents the volume density of a conserved 
quantity X and runs over the number density of a species, 
^a('":Oj the momentum density, p{r^t), and the energy 
density e{r,t), and where Jx{r,t) represents the corre- 
sponding fiux. This formulation is broad in that it allows 



for the presence of a term ax, which can account for "sink- 
source" mechanisms in the dynamics of conserved quanti- 
ties. For example, when X represent molecular numbers, 
ffx can be used to describe the kinetics of chemical reac- 
tions. In the case where an electric field E is applied on 
the system, which may contain molecules carrying charges 
{q^}, a-p = EY.aI^'^A and ffe = Y^j^q^fiAE ■ va may 
be used to model the effects of the electric field, where va 
denotes the velocity of species A. 

In principle, the law of angular-momentum conserva- 
tion should also be included. In standard hydrodynamic 
descriptions of conventional fluids, the canonical approx- 
imation is that each local fluid element has no intrinsic 
angular momentum. In the case of the type of membrane- 
fluid systems under our considerations, it would be ex- 
pected that physical situations exist where the approxi- 
mation is a valid one, and also that in other situations it 
no longer holds. But, we will work with the simpler cases 
where the approximation may be made. 

Based on the structure of the model system, x{r,t) is 
expressed as 

x(r,t) = x+(r,i)6'+(r,i) -f x"(r,t)e'"(r,i) 

+ [ dA{e,e,t)^{e,e,t)s{r-R{e,e,t)) , (4i) 

Jm 

where x^{r,t) represents the volume density of X in the 
"±"-bulk fiuid, M indicates that the surface integral is 
over the whole dividing surface, and finally, a;(^^,^^,t) is 
the surface density of the excess of quantity X. It fol- 
lows immediately from the concept of the model system 
and the assumption of local thermodynamic equilibrium 
that n^(r,t)'s, p'^{r,t), e'^^r^t) and s'^^r^t) satisfy the 
expressions of the bulk thermodynamics, Eq. p5|) and 
Eq. (|16|l . Thus, the two equations provide operational 
definitions of the corresponding chemical-potential fields, 
^^(r,t)'s, hydrodynamic velocity fields v'^{r,t), temper- 
ature fields, T'^{r,t), and pressure fields (r, t). Regard- 
ing the thermodynamic characterization of the dividing 
surface, a similar assumption is made: the non-equilibrium 
surface density of the excess energy, e(^^,^^,t), is still 
functionally related to the other surface densities, 
P('?^C^^), s(^i,^^i), and nA(^^^^^)'s as wefl as the 
shape field i?(^\^^,t) according to Eq. |(23J. The conju- 
gate variables defined thus by Eq. t;(^^,^^,f), 
'^(C^jC^j^)j and ^'^(C^, C^j *)'s are then considered as the 
velocity, the temperature and the chemical potentials of 
the dividing surface, respectively, and /rs(C^jC^ji) can be 
identified with the mechanical force exerted on the divid- 
ing surface. Consequently, Eq. and Eq. also hold. 

Similar to that of x{r,t), the model expression of the 
corresponding flux Jx{r,t) is given by 

Jx{r,t) = J^{r,t)9+{r,t) + J^(r, t)0- (r, t) 



dA 



M 



ji'He,e,t){n-V) 5{r-R) , (42) 
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id: 



where Jxir,t) represents the flux in the "±"-bulk fluid, be identified: 



It is one of the model statements that the functional de- 
pendence of Jxif,t) on the relevant bulk hydrodynamic 
state variables and bulk transport coefficients is known. 
Quantities C^, t) and j^^^ C^, t) account for the 

excess contributions due to the inhomogeneity in the real 
system. 

The presence of a non-zero C^, t) in the above 

model expression is only necessary when X represents lin- 
ear momentum. In that case jp^-* (^^, C^, t) can be identified 
as the surface flux of angular momentum. It includes the 
contributions from both the motion and the mechanical 
stress in the material, and is non-zero even when there is 
no motion. The reason for it is that a complete mechanical 
characterization of a membrane within the Gibbs model 
system requires both an effective surface-stress tensor, in- 
cluded in j^"-* (^^, C^, t), and an effective, non-zero bending- 

moment tensor, represented by a non-zero jl^\^^,£,'^,t), in 
order that there is a mechanical equivalence between the 
real and the model systems when their respective distri- 
butions of mechanical properties are integrated over the 
transverse dimension across the inhomogeneous region. 
A more detailed discussion of this issue will be given in 
Ref. 1201 • 

Inserting the model expressions, Eq. (|41(l and Eq. (|42() . 
in the conservation law, Eq. H4Q(I . carrying out the partial 
differentiations using Eq. (|39|l and performing some par- 
tial integrations lead to the following set of equations: 



;± 



d_ 

Dtx 



V J, 



3f - xc 



(jjf • n — x'^dtR ■ n) 
(j^ • n — x~ dtR ■ n) + (Jx 



(43) 



(44) 



dtR) -n + D^ (j« • , (45) 

(46) 



where a new differential operator with respect to t, 
Dt{») = ^/^/gdti\/g'), has been introduced. and x^ 
represent the boundary values of the bulk hydrodynamic 
quantities, i.e., the values of J {r,t) and x^{r,t) evalu- 
ated at the dividing surface, respectively, ax is the surface 
excess of the rate of generation/disappearance of quantity 
X associated with the sink-source term in Eq. (|40|l . 

Eq. H43|l is a reiteration of one of the model statements 
that fluids with the known bulk behaviours fill the regions 
on the two sides of the infinitely thin dividing surface. 
Eq. (|44|) represents the set of equations governing the dy- 
namics of the relevant excess surface quantities. Eq. H45|l 
and Eq. 146|l are simply conditions of self-consistency, im- 
plied in the Gibbs-model construction, on the normal com- 
ponents of j^*'-' and j'^x\ 

Based on Eq. H44|l . a tangential current within the di- 
viding surface, denoted by j", and two transverse currents 
entering/leaving the dividing surface, denoted by j^, can 



(0) 



Jx 

jx=''X' 

Eq. I|44(l can thus be written as 

Dtp = -Dajp + jp - jp , 
DtUA = -DajA + Ja - J'a + X! '^A^kC'^ , 



^ Jx ■ n — x^dtR ■ n 



(47) 
(48) 

(49) 
(50) 
(51) 



for cases where tJj, = Ue = 0. We will only consider such 
cases in what follows. 

The last term in Eq. (|51l) has been added to allow for 
the possibility of "chemical reactions" taking place in the 
membrane. The summation index K runs over all possible 
reactions, is the rate of reaction K per unit area and 
VA,K the stoichiometric coefficient of species A in reaction 
K . The word "chemical reaction" in our theory should be 
understood in a broader sense than that pertaining to a 
genuine chemical reaction. Connected to the assumption 
used in the formulation of thermodynamics that the mem- 
brane is considered to be quasi-statically impermeable to 
certain chemical species, non-equilibrium transport across 
the membrane of any of those chemical species (denoted 
by (7), such as the flip-flop process of a particular lipid 
species, is modelled by the following "reaction," 



(52) 



where C"*" and G~ are considered as two different labelled 
species. 

It must have not escaped the reader that Eq. H41|l 
and Eq. H42|l on their own are not sufficient to define 

x{e,e, t), j^°\e,e,t) and ji'\e,e. t). At the concep- 
tual level, the definitions of the surface excess quantities 
may be understood in the following sense. Consider the 
cell S defined in Section 3 and bear in mind in particu- 
lar that must be chosen such that the hydrodynamic 
behaviour of the real system coincides with that of the 
model system outside the top and bottom surfaces of the 
cell. If X£ denotes the amount of quantity X in the cell 
in the real system, the following condition of equivalence 
between the real and the model systems 



dV x{r,t) 



(53) 



then defines x(^^,^^,t). In a similar fashion, a number of 
conditions of equivalence should be satisfied by the model 
quantities jI"^ (^\ ^^j) and j^^^ (C\ t). If represents 
the cross section of S at constant then the following 
model quantity 



dA ■ [Jx - xdt {R + hn) 



(54) 



where dA is a normally directed area element on _B" and 
where the integration is taken over the transverse dimen- 
sion, ft,, along the normal n, should equal to the total 
current crossing B°' in the real system. 
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The discussion in the preceding paragraph in fact gives 
rise to a subtle issue concerning the relationship between 
the surface density of the excess momentum, p, defined 
by Eq. H53|l . and the excess current associated with the 
total mass, jp"-*. Given p, the hydrodynamic velocity v 
associated with the dividing surface is determined by the 
local thermodynamics, i.e.. 



V = p/p , 



(55) 



;(0) 



where p is the surface density of the excess mass, j 
defined by the condition Eq. (|54|l . is not equal to p in 
principle. The difference is on the order of /R where 
R is curvature, and it will be neglected, to a first approx- 
imation. In what follows, we will thus use 



4"' = pv 



(56) 



By this approximation, Eq. H45(l reduces to 

V ■ n = dtR ■ n , (57) 

a description consistent with our intuition. An alternative 
expression of Eq. lf57|) . which will also be used later, is 



V = v^ta + dtR , 



(58) 



where v°'ta can be interpreted as the surface-intrinsic part 
of V. 



5 Entropy production 

The set of equations of dynamics derived in the previous 
section, Eq. (gHl, Eq. and Eq. (OJ, involve both the 
currents describing the transport processes in the tangent 
space of the dividing surface, j", and the currents de- 
scribing the transverse transport processes, which in turn 
involve the boundary values of the bulk hydrodynamic 
fields. In order that the equations form a closed set, re- 
lations between the currents and the surface thermody- 
namic/hydrodynamic fields should be developed from the 
equation of entropy balance, according to one of the ba- 
sic principles of non-equilibrium thermodynamics. In this 
section, we present the derivation of the equation of en- 
tropy balance from which entropy production associated 
with transport processes is identified |24|. 

The tangential currents can be decomposed into two 
parts, a reactive part, which is associated with reversible 
transport processes, and a dissipative part, which is asso- 
ciated with irreversible processes: 



Ja ~ iA.r + 3aa 7 

• Q -a I -g 

J p J p,r ' J p,d 1 

■a _ -a I -a 
Je Je,r ~ Je,d ' 



(59) 
(60) 
(61) 



This is due to the fact that, when the dividing surface is 
curved, the volume measure and the area measure will then 
have a different dependence on the distance to the dividing 
surface. 



where the subscript "r" denotes the reactive part and "d" 
the dissipative. Once an expression for the entropy pro- 
duction is obtained, the reactive parts of the currents can 
be determined, by the definition that they should not con- 
tribute to the entropy production; moreover, and more im- 
portantly, thermodynamic/hydrodynamic "forces" driv- 
ing the dissipative currents can also be identified |25I26| . 

The equation of entropy balance in its general form 
can be written as 



(62) 



where j" represents the tangential current of entropy, jf^ 
are transverse currents, and is the density of entropy 
production. The derivation of an explicit expression of the 
equation starts from calculating, based on Eq. H24(l . the 
variations with time of all the thermodynamic quantities 
and isolating Dts as 



Dts 



1 r 

T 



Dte + -dtR-v- Dtp 



A 



fi^DtUA - (5" • dtR) 



(63) 



Further derivation can be carried out by using the con- 
servation laws, Eq. (|49|l to Eq. (|5HI . to replace the corre- 
sponding time derivatives and using the identity derived 
in Section 3.2, Eq. H38|) . to replace /^g • dtR with 



/rs ■ dtR= frs ■ - V^ta) 



/,s • ■» -f V" • daV + sdaT + ^ UAdaP,^ 



(64) 



where Eq. has been used. With the use of /^g = 

Z?c(S'"q) in addition, the equation of entropy balance can 
be rearranged into 



DtS = -D, 



}-(^j- + s'^-dtR-Y,p^rA 
r,+s^ -dtR-Y^ti^fA-^-jp 



3 s J s 



1 



(0) 



K 



^ {je - Jt) - ■ (jp - j + ) 

- Ja) ~ Us ~ it) 



T 



(65) 
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where Fk = f^^'^A.K is the affinity for reaction K. The 
interpretations of the different terms in the equation are 
made clear immediately by a comparison with the general 
form, Eq. Hfi2(l . The terms collected in the first pair of 
square brackets give the tangential current of entropy, 



rp \ Je 



S" ■ dtR - 



A 



(66) 



The terms in the second and third pairs of square brackets 
sum up the entropy production from all of the membrane- 
related transport processes. Those in the second pair rep- 
resent the entropy production from transport processes 
intrinsic to the dividing surface, and those in the third 
describe the entropy production from processes of trans- 
port between the bulk fluids and the dividing surface. 



5.1 The intrinsic transport processes 



the rest. Given that a judicious choice A = O can be made 
for one reason or another, Eq. (|72(l can be written as 



Jo,d 



A^O 



(73) 



The entropy production from the intrinsic dissipative 
processes can finally be expressed in terms of the various 
independent dissipative currents, j^^, Ja^o ^^'^ Jed 
and the thermodynamic forces driving them: 



jf,ddaT-j;^^-d^v 



A^O 



do 



K 



K 



where 



TO 



1 

V~9 



6E 
Sua 



p,p,s,R,{nB^A.o} 



(74) 



(75) 



The reactive parts of the currents can now be determined 
from the explicit expression of the entropy production 
given in Eq. (|65|l . An examination of the terms enclosed 
by the second pair of square brackets yields 

Jp,r = -^ro)+P«"' (67) 
jX, = UAv" , (68) 

and 




- 5" • dtR , (69) 



because these currents alone, in the absence of the dis- 
sipative parts, make no contribution to the entropy pro- 
duction. Finally, it can be concluded that the reactive and 
dissipative parts of jf defined in Eq.JBB)) are given, respec- 
tively, by 

Jf,r = , (70) 

j:,d = ^ (^3td - ^ ■ jp,d - E t'^^X^ ■ (71) 

In the most general sense, the above identifications of the 
reactive currents are not complete. We will discuss this 
point again where constitutive relations are derived. 

The various dissipative currents, j^^, j"^ ^'s, and ^ 
can now be determined. However, not all of the ^'s are 
independent due to a constraint 

E^'^^Id^O, (72) 

A 

which follows from Eq. ^ and Eq. Thus, j^^ of 

any species A can be chosen to be the one dependent on 



5.2 The processes of transport between the surface 
and the bulk fluids 

The processes of transport between the surface and the 
bulk fluids contribute to the total entropy production in 
the form of those terms contained in the last pair of square 
brackets in Eq. (|^ . It is clear that the contributions de- 
pend not only on the boundary behaviour of the bulk hy- 
drodynamic fields, but also on the surface hydrodynamic 
fields. A more illuminating expression of the contributions 
can be derived as follows. 

The values of the bulk currents evaluated at the divid- 
ing surface appear in Eq. H48|) which defines the transverse 
currents, j^, , j^'s and jf, and they are given, respec- 
tively, by 

J± = (e±+p±)7;±-T±.7;±+J±, (76) 

where is the viscous stress of the corresponding bulk 
fluid and Jg is the heat flux; 

4 =P^v±v±+p^\-T^ , (77) 
where I is the identity tensor; 

J± = n^v± + J± ^ ; (78) 

and 

Jtd = (J^ - f4jXd)/T± ■ (79) 
Inserting these explicit expressions into Eq. (|48|) yields 

^p^[n-{v±-v)v±]+p^n-n-T^ , (80) 
jf = P^v']. [n ■ {v± - v)] +p^n ■ V 

-n-J^-v±+Y,l^±3i + T±jt , (81) 
A 
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where n ■ v = n ■ dtR has been used. 

The sum of the entropy production from all the pro- 
cesses of transport between the surface and the bulk fluids 
can now be expressed as 



T(Ts,l_ = je -jt-V- {jp - jp) 
A 

A 

+3- (T_ - T) - j+ (T+ - T) 

1-2 1 - 2\ r ^ 



- 2^ + 2^ j'>T'-{v+-v) 



(82) 



In the above expression, the currents j^'s and n-{v± — v) 
are related by 



The newly-introduced quantities, and n±, are de- 

fined by 



/4 



— /;0 



A 

m n± 



(85) 



TO 



0± 



Ito"^ (.±-.)V(/±%iTO°"4) 



(86) 



Their physical interpretations become more evident when 
it is recalled that (^^ -I- ^m^v'^) and {^^ + ^m^v^) are, 
respectively, the bulk and surface chemical potentials of 
species A in the corresponding rest frames. 



6 Constitutive relations 

The expressions of the entropy production derived in the 
previous Section allow us to identify the conjugate "force" - 
current pairs associated with all the different dissipative 
processes. In this Section, we describe how physically mean- 
ingful relations between the currents and the forces are de- 
veloped, under the assumption that the non-equilibrium 
dynamics of the system is within the linear regime. To 
follow the standard terminology of statistical physics, we 
will call those relations constitutive relations in general 



^TO^j± =p±n.(t;±-t;) , (83) 

A 



6.1 Symmetry based classification of forces and 
currents 



which follows from J2a''^^'^a ~ P^'^^i Eq. (|48() and 
Eq. H57|l . Consequently, j^'s for all A values may be 
chosen as the independent currents, or alternatively, 
p'^n ■ {v± — v) and {j^,A ^ 0''=} may be used as the 
independent currents, where 0+ and 0~ denote the two 
judiciously chosen species, whose currents will be elimi- 
nated explicitly. 

Making the latter choice and rewriting Eq. (|82|l finally 
gives 



Tas,± = J- (T- - T) - jt {T+ - T) 

A^0± A^0± 

+ (— n • • n + 77_) n ■ (v^ — v) 
+ (n ■ ■ n ~ 77+) n ■ {v+ - v) 

+ {n-T+-Q[t^-{v+-v)]. (84) 



Eq. H74|) and Eq. (|84|l involve many different currents, such 
as ip,di Ja d'^' ^^'^■i many different driving forces, 
such as daT, daV, daP'^, etc..^ It is well known that a sin- 
gle current may be driven by several different forces. A 
systematic way to identify all the possible different forces 
driving a particular current is the following: first, to clas- 
sify all the forces and currents according to their behaviour 
under a group of orthogonal transformations, consisting of 
both rotations of the internal coordinate system and the 
inversion of the local normal vector to the dividing surface; 
and then, to determine whether the symmetry of the sys- 
tem allows for, or forbids, a force of one type, generically 
represented by F^, to drive a current of another type, de- 
noted Jj. Here, the Roman superscripts/subscripts label 
such classifications of forces and currents. 

The generic classes of behaviour of a quantity under 
the coordinate transformations consist of the following: 

^ In the formalistic sense what quantities are called forces 
and what are called currents are completely arbitrary. The 
convention we have adopted here conforms with either physical 
intuition or historical usages. 
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scalar, vector, and tensors of different ranks, which trans- 
form like a scalar, vector, and tensor with respect to in- 
ternal coordinate transformations, but which remain in- 
variant with respect to an inversion of n; pseudo-scalar, 
pseudo- vector, and pseudo-tensors, which are different from 
scalar, vector, and tensors only in that they change their 
sign under the inversion of n. 

For the intrinsic dissipative processes involved in 
Eq. (|74|l . the classification with respect to internal coordi- 
nate transformations is more essential. When applied this 
leads to the following conclusions: 

a) Genuine chemical reactions involve forces and cur- 
rents Jk, which are given, respectively, by 



F 



K 



-F, 



K 



A 



Jk^^. (87) 



Whether F^ and Jk are scalars or pseudoscalars de- 
pends on the precise nature of a reaction. However, 
when a chemical reaction process K* refers to the "flip- 
flop" from one side of the membrane to the other of 
molecules of a particular chemical species, as described 
by Eq. (|52|l . the corresponding force and current 



K' 



are clearly pseudoscalars. 
b) Heat conduction and diffusion in the surface involve 
forces, F^ and {F^,A ^ O}, and currents, jf^ and 
{Ja d' ^ 7^ *^}' which are 



F^ = -daT 
ja _ -a 



F^ = ^d^il 

ja _ -a 
•J A = JA,d 



r-.A 



(89) 



These forces and currents transform like vectors under 
the rotations of the internal coordinate system, 
c) The dissipative momentum transport involves forces 
and currents, whose representations need both a tensor 
and a vector, 

daV — {daV ■ tfj) + {daV ■ Tl) Tl 

= -{Fo^ptf^ + F(n),an) , 



P:{n)' 



(90) 



3p,d = [Jp,d ■ t j t/3 + [3p,d ■ ' 

= J ftp + r 

Thus, -^(71), a and Jp (^n) f'^^'iT^ the (pseudo) vector force- 
current pair. The tensorial part, Fa/3, is not a symmet- 
ric tensor in general and can be decomposed into three 
contributions, 



For the transport processes that give rise to the en- 
tropy production given in Eq. (|84|l . it is more meaning- 
ful to use linear combinations of the apparent forces and 
currents to generate new forces and currents that either 
remain invariant or change sign under an inversion of n. 
Before we discuss the new forces and currents, a change of 
notation will be made concerning our reference to species. 
So far, the definition of the labelled species, represented 
by A, is used to keep the derivation concise. In what fol- 
lows, the chemical identities of the molecular species will 
be explicitly referred to, in order that the physical inter- 
pretations of quantities associated with material transport 
become more apparent. Specifically, C will be used to rep- 
resent those molecular species, which have been assumed 
to reach chemical equilibrium across the membrane quasi- 
statically, while C will denote those molecular species 
which are assumed not to permeate the membrane quasi- 
statically. With this notation and with the choice that the 
"0" in the species expression 0* used in Eq. H84|) refers to 
one of the C species, the relevant terms in Eq. (|84() acquire 
an alternative expression 



(92) 



From Eq. H84|l the new forces and currents can now 
be derived for the various different transport processes, as 
follows. 

a) The first four lines of the equation sum up the con- 
tributions from the absorption and conduction of heat 
by the surface as well as from the adsorption onto, 
and permeation across, the surface of the constituent 
molecules. Genuine scalar forces F(s-) and currents J^^\ 



as well as pseudoscalar forces F(a 
are given by 



and currents J^^\ 



T- - 2T) 



(93) 



(94) 



F, 



af3 



+ Eape'^^F^s , 



9o.pF\) 



(91) 



each of which is invariant under any internal-coordinate 
transformation and each of which should be considered 
as an independent force. The corresponding currents in 
Jap can be identified in a similar way. 



f: 



1 

2 



c 



(a) 



2 (-^^ 



-Jc) 



(95) 



(96) 
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c ^ 2 ^ c+ 



c- 



1 ..^ 



2 ^-^CH 



7-(a) 



7(a) _ 



1 



2 ["•(T^ + Td)•"-/7+-ff- 
7^ • («+ - v^) , 
^["•(Tj-TZ)-n-77, 
n • (t>+ + — 2t>) . 



(97) 
(98) 

(99) 

) 

(100) 



The interpretations of the scalar currents and the pseu- 
doscalar currents are rather obvious: The scalar cur- 
rents describe the heat absorbed, the material adsorbed, 
and the pseudoscalar ones describe heat conduction, 
material permeation across the surface, 
b) The last two lines of Eq. H84|l describe the dissipation 
associated with the hydrodynamic "slip" between the 
surface and the two contacting bulk fluids. The forces 
and currents involved can be written as genuine vectors 
and pseudovectors, 



(101) 



(S),Q 


— 




2v) 


7(s),a 

Up 


1 

— —n 
2 




) -t" 


Ff. 

(a),Q 


— 




7 


7(a), a 
dp 


1 

— —n 
2 




) -r 



(102) 



6.2 Constitutive relations 



The total entropy production associated with the dividing 
surface, the sum of Eq. H74|) and Eq. H84|l , can now be ex- 
pressed in terms of the forces and currents defined above: 



AF' 



(103) 



where i runs over all the classes listed above. The basis for 
developing constitutive relations is our requirement that 
Tas be positive definite, although it is developed for the 
dividing surface in the Gibbs model system, which is not a 
physical surface. This requirement ensures the thermody- 
namic stability of any equilibrium state described by the 
theory; its validity will be discussed in the last section of 
the paper. 

The most general form of the constitutive relations is 
given by 



Ji — ^ ^ ^ijF'^ , 



(104) 



where J7y are phcnomcnological parameters. It immedi- 
ately follows from Onsager's reciprocal relations that 's 



must satisfy a general property: i7.y = Qji if Ji and Jj 



have the same parity under time reversal, and 



if Ji and Jj have the opposite parity. The positive defini- 
tiveness of Tos is then ensured if matrix flij is positive 
definite. The anti-symmetric elements of flij require a few 
more words. It is not difficult to see that they do not con- 
tribute to the entropy production, in other words, any pos- 
sible couplings characterized by anti-symmetric J7y are in 
fact reactive couplings rather than dissipative couplings. 

In principle, any of the currents, Ji, may be driven 
by any of the forces, Fj, leading to the so-called cross- 
coupling. In other words, in the absence of pertinent sym- 
metries or invariances, all forms of cross coupling are pos- 
sible. However, if the physical description of the system 
is invariant with respect to some or all of the orthogonal 
transformations, then the invariance will eliminate certain 
cross-couplings. For example, if the equilibrium shape of 
the membrane (or the dividing surface) is symmetric with 
respect to rotations around any local normal vector, then 
scalar, vector and tensor forces can only drive currents of 
the same type, leaving the coefficients of all the related 
cross-coupling terms zero. Similarly, if the physics of the 
system does not distinguish one side of the membrane from 
the other, a pseudoscalar force can not drive a scalar cur- 
rent and a scalar force can not a pseudoscalar current. 
Consequently, the corresponding J7ij's vanish. If it turns 
out that these symmetries in practice are not exact, but al- 
most correct, then the corresponding cross-couplings will 
be weak. An obvious approximation is to discard those 
cross-coupling terms. It may be worth noting, however, 
that symmetry properties alone are not sufficient for iden- 
tifying physically meaningful and relevant forces, currents, 
and constitutive relations. The general physics embod- 
ied in the entropy-production equations and the specific 
physics of a system under consideration are necessarily 
needed. 

A few explicitly written constitutive relations may lead 
to some more intuitive understanding of the general dis- 
cussion given above. For example, concerning the surface- 
intrinsic processes described by Eq. (|89|l . there are 

J," = ifT(-9"T) + a^^(-9"A^) , 

Jl=Y. ^?A,s(-a"/i^) + QaA-Q'^t) . (105) 

B^O 

These two relations should be familiar to the reader, where 
Kt may be interpreted as the effective heat conductivity 
of the surface, and /2a,b's are the effective dissipation co- 
efficients associated with diffusion. Another example can 
be developed based on Eq. (|90l) and Eq. (|91l) , and has the 
following form: is the following: 

- C {dpv ■ t^) r , (106) 

where the phcnomcnological constants Q and ?7s may be 
interpreted as the only relevant surface viscosity coeffi- 
cients. This constitutive relation is a two-dimensional ana- 
logue of the familiar expression of viscous stress in a bulk 
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fluid. To be sure, the symmetry reasoning alone would al- 
low for the presence in the above constitutive relation of 
both the last antisymmetric part of the force in Eq. H91|l 
and a "force" term of the form {d°'v ■ n) n. The physi- 
cal requirement that no dissipation should be associated 
with uniform rotational motion of the whole system, how- 
ever, renders those two terms absent. The latter term was 
allowed to be present in the work reported in Ref. |30| . 
though. 

The processes of transport of material across the sur- 
face are also of practical interest. One of such processes is 
the so-called flip-flop process. In the real world of mem- 
brane systems it often occurs that this process is so slow as 
to be practically irrelevant for experimental observations 
|18| . However, situations where this process is relevant are 
also encountered j3lj . where the kinetic rate characteriz- 
ing the process can be determined experimentally. In the 
framework of our theory, the simplest form of the consti- 
tutive relation associated with the flip-flop of a chemical 
species C is given by 



K' 



(107) 



where the phenomenological parameter fl^ may be related 
to experimentally determined flip-flop rate. A related con- 
stitutive relation concerns the permeation of molecules 
of a C type. As it has been pointed out in connection 
with Eq. 1)98(1 . the current of permeation, J'q \ is a pseu- 
doscalar, thus can be driven by F'j\^ as well as by = 
— {pP^ — fjp ) defined in Eq. 1)88(1 . The corresponding con- 



stitutive relation may be written as 



C,pcrm 
C -pcrm 



(108) 



In contrast, in the constitutive relation characterizing the 
permeation of molecules of a C-type, only the counterpart 
of the dissipative coefficient f?'^' exists. 

^ C,perm 

The final example concerns the total permeation of 
material across the surface described by the forces and 
currents defined in Eq. ((100(1 . If it may be assumed that 
there is no cross coupling between the scalar and the pseu- 
doscalar quantities, then the following simple constitutive 
relation may be written: 

n ■ {v+ + w_ — 2v) 

= /?perm^ [n ■ (T+ - T^) • n - 77+ + 7T_] . (109) 

The constant J7perm is related to membrane permeability. 
This relation would be the familiar Darcy's law [32, but 
for the presence of the viscous stresses from the bulk fluids. 



6.3 From constitutive relations to equations of motion 

In the above, we have presented a framework of principles 
for developing constitutive relations, the specific forms of 



which must depend on the particular physical system un- 
der consideration. We end this whole section with some 
comments on how the general theory may be considered 
from an operational point of view. 

The complete set of equations of motion may be di- 
vided into two subsets: those that describe the hydrody- 
namics of the bulk "filler" fluids, i.e., Eq. ((43(1 . and those 
that describe the dynamics of the surface excess fields, 
i.e., Eq. (O. Or, more specifically, Eq. (gU to Eq. ifCT) . 
together with the one that governs the conformational dy- 
namics, Eq. lf57|l . or Eq. lfS5|l . These two subsets of equa- 
tions can only be solved if specific constitutive relations 
are given or developed. The constitutive relations associ- 
ated with transverse transport processes provide the nec- 
essary coupling between the boundary values of the bulk 
hydrodynamic fields and the surface dynamic fields, as 
Eq. ((109(1 illustrates, for example. They may be viewed 
and used as boundary conditions for the bulk equations of 
motion. An operational strategy may then be to solve first 
the bulk equations of motion and express all the boundary 
quantities of the bulk fiuids appearing in the surface equa- 
tions of motion in terms of the surface dynamic fields, and 
then obtain an effective, closed set of equations of motion 
for the surface dynamic fields. As the focus of a study of 
dynamics of a membrane is naturally on the surface fields, 
the effective equations should be the basis for analyzing 
the dynamics of the membrane. 

The number of the final, effective equations of mo- 
tion, or the number of relevant dynamic fields, associated 
with the dividing surface can be worked out in a gen- 
eral fashion. Similar to the case of a bulk fiuid, relevant 
conserved variables are given by the surface densities of 
excess energy, excess momentum, and numbers of excess 
particles, represented by e, p and {^a}! respectively, as in 
Eq. 1(49(1 to Eq. 1(51(1. An extra, non-conserved, and scalar, 
dynamic variable, however, becomes necessary for describ- 
ing the conformational aspect of membrane dynamics, as 
mentioned in the preceding paragraph ^. For a reader with 
a particular interest in the linear analysis of equations of 
dynamics, the total number of the relevant dynamic fields 
in principle determines the number of independent dy- 
namic modes that are implied by the effective equations 
of motion. In practice, the number and the specific na- 
ture of those modes that are relevant to, or observable in, 
experimental observations of the dynamics of a particular 
system will have to be determined on a case-dependent 
basis. 

Given its importance and accessibility for experimental 
observations [2| the conformational aspect of membrane 
dynamics on its own merits a comment. Eq. ((57() may 
be considered the corresponding equation of motion, and 
clearly, it is coupled to the other equations of motion. 
That coupling is expected to be complicated in a general 
case. In some limit cases, however, an explicit form of the 



^ In principle, the shape field, ii is a vector, implying three 
scalar variables. But, only one variable, corresponding to the 
motion of the membrane in its normal direction, is the phys- 
ically relevant one, as a consequence of the reparametrization 
invariance [8()( . 
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equation is available. One of such limit cases, the case of 
membranes made of single-component lipid-bilayers, has 
been explored earlier 7; and will also be discussed later 
in Sec. 17.21 of this paper. An effective equation for the 
conformational dynamics in that limit case can be found 
in|7|. 

7 The general theory in a few limit cases 

The formulation of the theory as presented in the previous 
sections is kept very general. This feature is useful, seen 
from the formalistic point of view. A few limit cases are, 
however, interesting from the practical point of view. In 
fact, another prediction of our theory, which is in principle 
non-trivial, will be brought out more clearly in one of the 
limit cases. In this section we will discuss these limit cases. 

7.1 Multi-component membrane with no-slip boundary 
conditions 

In treatment of the bulk hydrodynamics of a viscous fluid 
in contact with a boundary, the boundary condition that 
is usually imposed is the so-called no-slip boundary condi- 
tion, where the bulk fluid at the boundary is constrained 
to have the velocity of the boundary material. 

If this limit may be assumed as a valid one for the 
membrane-fluid system under our consideration, its im- 
plementation does not appear straightforward, due to the 
peculiar structural characteristics of the membrane. On 
the one hand, one type of molecular building blocks of 
the membrane are amphiphilic molecules, which are dis- 
tributed in two opposing monolayers, each in contact with 
a different bulk fluid. There is no physical reason to as- 
sume that the lipid molecules forming the two different 
monolayers should have the same collective velocities, and 
in turn, no reason to assume that the boundary velocities 
of the two bulk fluids are the same. On the other hand, 
another type of the molecular building blocks are trans- 
membrane proteins, which span the whole thickness of the 
membrane and which contribute to the dynamic interac- 
tion between the membrane and both of the bulk fluids. 
The question is, in other words, what surface velocities 
should be matched tangentially with the boundary veloc- 
ities of the two contacting bulk fluids. The previously de- 
fined surface velocity v, which is the center-of-mass veloc- 
ity of the whole surface, obviously is not the appropriate 
one. 

We propose, as an answer to the question, that the no- 
slip boundary conditions be implemented by matching the 
boundary velocities of the two bulk fluids with two in prin- 
ciple different linear combinations of the velocities of all 
the different species associated with the dividing surface. 
By use of the fact that the collective velocity of a partic- 
ular species A is given by haVa = Ja = (j'S^d*" + nAvj , 
the no-slip boundary conditions can be formulated as 

V±^Y.^±OAAto. + nAv) . (110) 
A 



The L^'s in the expression are phenomenological quan- 
tities, and may be interpreted as the effective areas of 
the different species "seen" by the bulk fluids. Note that 
-^±"-^ = 1 must be satisfied in order that the expres- 
sion also hold for the case of uniform motion of all species. 

Using this constraint together with Eq. H72|l leads to an 
alternative expression of the no-slip boundary conditions 

v±=v+J2 LijAA^c , (111) 

A^O 

where 

Li^Li- . (112) 

Eq. (|lll|l implies that • {v+ - v) and t" ■ {v- - v) 
are no longer independent of the currents {Ja 0}-> 
and it can then be used to eliminate t" • {v± — v) from 
those terms in Eq. (|84|l that contain them. Collecting all 
the terms in the total entropy production which contain 
{j2 d' ^ 7^ including those in Eq. ((7^ . finally leads to 
the identification of 

= -do.fi'' + • T+ • t„ - L^n ■ • t„) (113) 

as the thermodynamic force conjugate to the diffusion cur- 
rent J2- 

The terms contained in parenthesis in the above equa- 
tion represent contributions from the two bulk fluids, and 
they underscore one of the specific predictions of our the- 
ory: the exchange of momentum between the surface and 
the bulk fluids is coupled to the dissipative processes in 
the surface. At the level of principles, this is a non-trivial 
prediction, and is made, to our knowledge, for the flrst 
time here. In Section |S1 an order-of-magnitude estimate 
of this new effect will be made with a view towards its 
experimental veriflcation. 



7.2 Single-component lipid bilayers 

The case where the membrane is a bilayer composed of a 
single species of lipid molecules is also worth discussing 
in connection with the earlier works dealing with such 
a system In those earlier works, the description of 

the hydrodynamics of the membrane was established in 
an intuition-based manner rather than from a systematic 
derivation. In this section, we will demonstrate that the 
earlier description can be formulated as a limit case of 
the formalism developed here, thus putting the earlier de- 
scription on a flrmer basis. 

In the earlier description, a number of limiting assump- 
tions were made. First, the temperature was assumed to be 
uniform in the whole system; secondly, the membrane was 
assumed to be impermeable to solution (water) molecules; 
thirdly, it was assumed that no transport of the lipid 
molecules took place either between the bulk fluids and the 
membrane or between the two distinct monolayers of the 
membrane; fourthly, it was assumed that the intramono- 
layer shear dissipations could be neglected, but that the 
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dissipation associated with the relative motion between 
the two monolayers was relevant; finally, the Stokes con- 
dition was assumed, i.e. that the inertia of the membrane 
surface could also be neglected. Consistent with these as- 
sumptions, the single-component lipid bilayer was consid- 
ered as a composite of two systems (+/—), one for each 
monolayer. Consequently, the description of the dynam- 
ics became much simplified, consisting of two continuity 
equations 

Dtni^ + D^j^^ = Dtm^ + {ni^W^) = , (114) 

where vi^ = W^ta + dtR were the velocities of the lipid 
molecules in the two monolayers, and by two additional 
equations which correspond to the conservations of the 
monolayer momenta under the Stokes condition, 



J n 







(115) 



/j-s represented the restoring force acting on each (+/— ) 
monolayer, which should be derived from the thermody- 
namic free energy, and was the force acting on one 
monolayer by the other. As a model for the intermono- 
layer dissipation, the following phenomenological expres- 
sion was employed 

• ft. = -bt" ■ {vi^ -vi^)= -b (W^ - W^) , (116) 

where b was the "intermonolayer friction coefficient." In 
order that the equations of the dynamics be completely 
closed, intuition-based expressions were proposed for f^,, 
which satisfied 



and 



J IS ''o 



J r 



(117) 



(118) 



f rs f rs • 

Eq. (|115|) can now be expressed in an alternative form, 
which will make a comparison with the current formalism 
easier. The sum of the two subequations yields the equa- 
tion corresponding to the momentum conservation for the 
whole bilayer, which is given by 



T^ 



T- 







(119) 



and the subtraction of the two subequations leads to the 
following expression 



Ji_.d 



?^/+ + ni_ 



1 

ni 



-1^ 



1 

ni_ 

a,.l-\ 



-T 



where 



Jl+,d = Ju 



ni± 



ni+ + ni_ 



(120) 



(121) 



obviously satisfying j^^^^ = -ii_,d- 

When the same system is considered within the cur- 
rent formalism under the same limiting conditions, the 
equations of dynamics for the system have the same for- 
mal expressions as Eq. (I114II and Eq. (|119|l . where jf^ ^ 
defined by Eq. H121|) represent the "diffusion currents." 
If, in addition to all the limiting conditions already men- 
tioned, the limiting, no-slip boundary conditions Eq. I|lll|l 
are also used, which in this simple case reduce to 

v+ = vi^ , V- = vi_ , (122) 

or equivalently to 

4+ = 4+ = ^ , =-l'- =-^, (123) 

the constitutive relation concerning the "diffusion cur- 
rents" is given by, following Eq. (|113|l . 



i7r 



— T+ - ^T- 1 . r 

ni_ 



5a(M'+-/i'-) 



(124) 



It is obvious that, if the following identification is made. 



ni+ + ni_ I b 



(125) 



the above constitutive relation becomes equivalent to 
Eq. (|12()|l . The above equation makes it clear that the term 
"diffusion current" within our formalism has a broader 
meaning than its canonical one, in that it describes rel- 
ative motion between species and does not have to be 
associated with processes of molecular mixing. 



7.3 Incompressible membrane and bulk fluids 

Theoretical work dealing with a generic membrane system 
often in practice employs the assumption that both the 
bulk fluids surrounding the membrane and the membrane 
itself are incompressible. This assumption is justified in 
most cases. With a view of applications of our general 
theory to those practical situations, we in this subsection 
discuss the limit case of volumewise incompressible fluids 
and membrane. 

In what follows, we will model the limiting condition 
of incompressibility by assuming that the average molec- 
ular volume associated with each molecular species in the 
membrane-fluid system remains constant throughout the 
system, in other words, being insensitive to the spatial 
inhomogeneities in other thermodynamic variables. In ad- 
dition we will only consider the case where another lim- 
iting condition holds, that there is no transport of any 
molecules across the membrane. This latter condition im- 
plies, according to the definition of "species" established 
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earlier, that in the Gibbs-model description of the system, 
two labelled species, (7+ and C~ , are associated with a sin- 
gle chemical species C. It is not difficult to see that these 
conditions amount to two separate constraints on the sur- 
face fields {ng± }. For our purpose, the two constraints are 
expressed in the following form 



«o± = B± {{nc+^a+,nc-^o-}) 



(126) 



where C = represents the molecular species (water) 
forming the solvent. 

The immediate consequence of the constraints is that, 
among all the transverse currents of material transport, 

(^{j^±^Q±}i P^'^T' ■ {'"± ~ ■")) : two - one for the -region 
and one for the "—"-region - must become dependent 
on the rest as well as on the surface currents, {jg± j}- 
There are different ways to remove these dependent cur- 
rents from the entropy-production equation. An obvious 
one is by direct substitutions of their explicit dependences 
on the other currents. Although straightforward, the im- 
plementation is tedious and changes the formulistic struc- 
ture of the derivations we have described so far. As a for- 
muHstically much more concise and simpler alternative, 
we have developed a different approach to implement the 
constraints, which is sketched below. 

In our approach, the following energy-density function 
is introduced first 



e(no+ , no- , {"c+5^o+ > '^c-^o- }) 

= ephys({n(5+-^o+''^C-#0-}) 

+ A*'^ (no+ - B+) + (no- - B_) 



(127) 



where ephys({n(5+-^o+ j '^c-^^o- }) represents the physical 
energy- density function for the constrained system, and 
where A*^ are two multipliers to be determined. It is 
obvious that, when the constraints are used, e(no+,no-, 
{n(j.+ /o+ , n(j.-_^o- })i coincides with the physical function 
ephys({n(5+^o+: '^(7+5^0+}) no matter what values the two 
multipliers take, and this fact will be taken advantage of 
in what follows. 

In the derivation of the entropy production the sys- 
tem is then treated formally as if it were not constrained 
by letting e(no+ , no- , {nc+-io+ > "c+^o+ }) pl^y the role of 
the surrogate of the physical energy-density function. The 
final expression of the entropy production. 



Ta. = 



r.A 



j7(r_-r)-j-+(r+-T) 



K 



+ (— n • • n + i7_) n ■ {v^ — v) 
+ (n • Tj[ • n - 77+) n ■ {v+ - v) 



(128) 



thus looks formulistically identical to the sum of Eq. (|74|1 
and Eq. (jHJ), except that the multipliers A" appear in 
the places of the surface chemical potentials /x" in 



m 



(129) 



n± r 



774 



(A"* + irn%2) 



(130) 



Eq. H128|) still contains all currents, both independent 
ones and the corresponding dependent ones. The freedom 
in choosing the values of the multipliers A'^ when the 
constraints are satisfied now means that, given a particu- 
lar choice of the independent currents, the corresponding 
dependent ones can then be eliminated from Eq. 1)128(1 by 
choosing values for A° accordingly. For example, for the 
choice where ({j^± , C ^ 0}) are used as the independent 

transverse currents, the values of A° must be chosen such 
that 



n • TJ • n - 77+ = 



(131) 



This choice renders both the fourth and the fifth lines in 
Eq. 1(128(1 zero, thus eliminating the corresponding depen- 
dent currents p^n ■ {v± — v). 

The values of A*^* chosen as such, together with the 
constitutive relations which can be derived consequently, 
should then be used in conjunction with equations of mo- 
tion where the constraints are explicitly satisfied, in order 
that the whole formulation of the dynamics form a self- 
consistent and closed one. 

The implications of Eq. ((131(1 merit a short comment. 
It follows from the the definition of 77± that the values 
of A° are linearly related to the chemical potentials of 
the bulk solvent molecules, which in turn are related to 
the pressures in the bulk fluids. Thus, bulk-pressure gra- 
dients tangential to the surface, ta ■ Vp± , are translated by 
Eq. I(131() into tangential gradients of A° and can thereby 
drive diffusion currents within the surface in principle, as 
implied by Eq. ((113(1 . The practical relevance of this mech- 
anism can, however, be argued to be insignificant. A quick 
estimate reveals that the tangential gradients of A° are 
given roughly by w'^ta ■ Vp± , where w° denotes the molec- 
ular volume of the solvent. When compared with the con- 
tributions in Eq. ((1131) from the bulk shear stresses, this 
effect is suppressed likely by the ratio between a molec- 
ular length and the wavelength characterizing the bulk- 
pressure gradients, thus probably insignificant. This sup- 
pression means that, although the excess quantity of the 
solvent molecules associated with the surface may not be 
negligible, its effect on the motion of other surface-related 
molecular species may be assumed to be insignificant. Our 
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understanding here justifies, therefore, the use of this as- 
sumption in the previous works on membrane hydrody- 
namics ,7.10) . 

A general remark may also be made to conclude this 
subsection. Although our method of constraint implemen- 
tation has been applied in the specific limit case of sys- 
tems consisting of incompressible fluids and incompress- 
ible membranes, it may also be applied to cases where 
constraints are different. The difference in the use of the 
method will only result from the different details of the 
constraint functions used in Eq. H126|) . 



8 Discussions 



In the previous sections, a formulation of hydrodynamics 
of the membrane-fluid systems has been presented in a 
very general form, both in terms of equations of motion 
based on the basic conservation laws and in terms of de- 
rived constitutive relations. Although an application of the 
general theory to any specific system is beyond the scope 
of this paper, we will, nevertheless, discuss in this sec- 
tion a number of general issues which must be dealt with 
in applications and which will also make the connection 
between the theory and experimental situations visible. 
For the following discussion we reiterate that the theory 
is only an effective one and that its purpose is to give 
a mesoscopic/macroscopic-scale description of the mem- 
brane dynamics, which can be connected to experimental 
descriptions on the corresponding scales. 

The first issue concerns the structure and the content 
of the theory. From the way the theory is constructed, it 
is clear that its description of the motion in the regions 
close to the membrane interface is, by design, not correct 
qualitatively. The obvious questions are, then, how the 
quality of the theory can be controlled, and what crite- 
ria should be used as measures of the quality of the the- 
ory. The means of the control is already contained in the 
theory itself and rests on tuning the set of kinetic coef- 
ficients associated with the various constitutive relations. 
Whenever the theory should be applied, it is assumed that 
the hydrodynamic conditions at the boundaries of the en- 
tire system are controlled by macroscopic means and are, 
therefore, known. Given those boundary conditions, and 
provided that the dynamical behaviour of the system in 
the regions sufficiently far away from the membrane (or 
the dividing surface), or the "bulk behaviour," is known 
already, or may be known, by "measurements," the set 
of kinetic coefficients can then be tuned in the process of 
solving the equations of motion under the given boundary 
conditions, such that the quantitative description of the 
"bulk behaviour" given by the theory matches that given 
by the measurements to some desired degree of agreement. 
Moreover, in the theory, a few quantities pertaining to the 
interfacial region can be calculated. 



where S is the familiar volume element already used in 
Section 18. II and 

-e^ <h<e+ , = ^" ± zi^"/2| (132) 

represent the side surfaces of S and dA is the area element 
on this surface times its outward-pointing normal vector. 
The counterparts of these quantities, or at least some of 
them, in the real physical system may also be obtained 
from measurements. One example would be the measure- 
ments of distribution and motion of membrane proteins 
within the membrane by the use of modern techniques 
such as Fluorescence Confocal Microscopy and Fluo- 
rescence Correlation Spectroscopy |34) . Such experimen- 
tal information on the "surface behaviour" of the system, 
whenever available, must also be used to tune the theo- 
retical kinetic coefficients, thereby to control the quality 
of the theory. 

Besides tuning the kinetic coefficients such that the 
quantities labelled 1. and 2. above are identical for the 
Gibbs model and the experiment, the membrane param- 
eters should also be tuned such that the position of the 
shape R matches the position of the dividing surface in 
the experiment. 

The second issue is a general one concerning the tun- 
ing of the kinetic coefficients. It has been a requirement, 
made in Section 6.2, that the entropy production in the 
theory be positive definite in order that the equilibrium 
states of the theory be stable. One consequence of this 
requirement, among others, is that kinetic coefficient f^u 
for each i is positive. It turns out that in certain situa- 
tions, this stability requirement can only be fulfilled by 
an approximation in the tuning of the kinetic coefficients. 
The following example of permeation of molecules across 
a membrane illustrates the point. 

Consider the following simplified picture of a real (in 
contrast to the Gibbs model) membrane-fluid system. The 
membrane is taken to be a homogeneous slab of material 
of thickness 2d. A certain molecular species has different 
diffusion constants, D' and D, in the membrane region and 
in the two bulk fluids, respectively. At a distance Li from 
the membrane slab, there is a reservoir of the molecules 
with a constant chemical potential /ii, and similarly on 
the other side, a reservoir of chemical potential /i2 < /ii 
is placed at a distance L2 from the membrane slab. 

The quantity of interest is the steady-state molecular 
flux, j"'^'^^, between the two reservoirs. It is easy to work 
out that 



J 



real 



DD' 



D'{Li+L2) + 2dD 



(a*i - M2) 



(133) 



1. f^dV x{r,t) , 

2. /g„ dA ■ [Jx - xdt (R + hn) 



In the theory, the real system is replaced by a Gibbs 
model system: two bulk fluids, where the molecular species 
has diffusion constant D, separated by an inflnitely thin 
dividing interface with an effective kinetic coefficient f2 for 
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permeation of the molecules, which must be tuned. The 
two reservoirs of chemical potentials /ii and ^2 are then 
at distances Li + d and L2 + d from the interface. In order 
that the steady-state molecular flux be reproduced in the 
model system, it is necessary that the boundary chemical 
potentials of the bulk fluids at the interface be different, 
denoted by fj,~^ and /i~ . The steady-state flux, j 
then be expressed as 



model 



Ddcl 



= -D 



-D 



- ^J■^ 

L2 + d 

^1- - fil 



(134) 



Li+d 

From this equation an alternative expression of can 
be obtained 

on 



J 



model 



n{Li +L2 + 2d) + D 
In order that the requirement, j™°doi 



(M1-M2). (135) 



filled, i7 must be set to be 



DD' 



2diD-Dr ^'''^ 

Eq. 1136() immediately illustrates a problem of the the- 
ory in the case where the molecules diffuse faster than they 
do in the bulk fluids, i.e. D' > D: on the one hand, the 
effective kinetic coefficient fi must be given a negative 
value; on the other hand, any equilibrium state becomes 
unstable in the theory when Q is negative. 

This problem does not, however, render the theory in- 
applicable to the case. A reasonable approximation can 
still be made, which involves setting fi to infinity. The 
chemical boundary condition for the bulk fluids at the 
dividing surface thus becomes ji^ = /i^. The error intro- 
duced in this approximation has an upper bound, which 
can be easily calculated: 

D' -D 



J 



real 



J 

real 



model 



< 



2d 
Li +L2 

2d 
L1 + L2 



2d 



D' 



(137) 



Clearly, the theory should still give quantitatively reason- 
able result if d is much smaller than Li or L2, or in a more 
general situation, if the length scale associated with the 
chemical gradient in the system is much larger than the 
microscopic thickness of the membrane ^ . 

The negative permeability is a direct consequence of invok- 
ing the Gibbs model in cases where D' > D. An alternative 
formulation may also be envisioned, as suggested by one of the 
referees of this manuscript. In this formulation, the thickness of 
the membrane would be taken to zero without corresponding 
"compensation" in the theory for the effects that are associated 
with the finite thickness of the membrane. As a consequence, 
the kinetic coefflcient associated with transmembrane perme- 
ation stays positive in all cases. This zero-thickness treatment 
is obviously an approximation, as is our suggested solution to 
the problem of negative permeability. The errors introduced 
due to this approximation may be comparable to those associ- 
ated with our formulation. 



The last issue to be touched upon here is the ques- 
tion of how the various kinetic coefficients in the effec- 
tive theory may be related to experimental data on the 
corresponding transport processes. The precise answer to 
the question depends on the nature of the experiments 
that are done. For example, if a permeation experiment 
is done in the way described above, where the molecular 
fiux is measured, then an effective kinetic coefficient 
can be derived based on Eq. (|135|l . Unfortunately, very 
few experiments have been done so far on the investiga- 
tions of membrane-involved transport processes; concrete 
situations where this issue can be addressed specifically 
are hard to find. We hope, however, that our theory may 
provide some guidelines for designing useful experiments 
regarding membrane-involved transport processes and for 
interpreting data obtained. 

It would add strength to a general theory such as 
the one presented in this paper if it makes conceptually 
non-trivial and experimentally relevant predictions of new 
physical effects. Indeed, our theory does make predictions 
that are conceptually non-trivial, among which the fol- 
lowing constitutive relation - the simplest one that can 
be written down based on Eq. (|113|l . 



(138) 



is the most notable. An order-of-magnitude estimate of 
the effect predicted by Eq. H138() . namely, that the vis- 
cous stresses exerted by the bulk fluids on the membrane 
can drive the surface diffusion processes in the membrane, 
may shed some light on the experimental relevance of this 
theoretical prediction. 

Consider the case where a membrane is composed of a 
lipid species and a protein species embedded in the lipid 
bilayer. The subject of our interest is the motion of the em- 
bedded protein. Assume that the physical conditions are 
such that there exists a surface gradient of the protein con- 
centration in the membrane, sufficient to drive appreciable 
protein diffusion. In addition, the membrane is subjected 
to a shear flow induced in one or both of the bulk fluids. 
For an order-of-magnitude estimate it suffices that we as- 
sume near-plane geometry for the membrane and that we 
include in the chemical potential of the protein only the 
contribution from entropy of mixing, i.e. 



//P = kTlnrip 
It follows immediately then that 

kT An-r 



Up Al ' 



(139) 



(140) 



where Al denotes the wavelength characterizing the con- 
centration gradient. The contribution from the bulk shear 
stress can also be estimated, based on Eq. (|138|l . to be 



(141) 
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where is in order of magnitude similar to defined by 
Eq. (fTn|l with the species index O referring to the lipid, 
but in explicit expression different from Lj_ (19) . 

For the newly predicted effect to be comparable to the 
effect associated with the chemical gradient, the variation 
of the velocity characterizing the bulk shear flow should 
be 

Av^^^. (142) 
r]LP Up 

Taking 7/0 = lO^'^Pa • s, the shear viscosity of water at 
room temperature {kT = 4 • 10~^^J) and setting LP ^ 
500A^ = 5 • 10~^^m^, a reasonable value for the cross- 
sectional area of a typical transmembrane-protein mole- 
cule, we get 



Av . 



Atlr, 



Vo 



(143) 



Based on the above expression, it may be concluded that 
it is very plausible that protein diffusion induced by bulk 
shear flows can be observed experimentally, if the shear 
viscosity (or viscosities) of the bulk fluids can be increased. 

Finally, we conclude this paper with our outlook, from 
the stand point of this work, on study of non-equilibrium 
dynamics of membrane systems. It is our belief that non- 
equilibrium behaviour of membrane systems must be in- 
vestigated, described and understood. It is then our hope 
that the theory presented in this paper will contribute 
to that study as a general framework, from which applica- 
tions to specific concrete problems can grow. Compared to 
more intuition-based approaches to formulating descrip- 
tions of non-equilibrium dynamics of membrane systems, 
the present formal approach has the virtue that it avoids 
guesswork and ensures that relevant effects or mechanisms 
will be included and correctly described in the equations 
of motion and in the boundary conditions. To be sure, 
not all of the constitutive relations rendered possible in 
the general theory are relevant in a specific application to 
a specific system. But, a systematic analysis of what is rel- 
evant or what is not can be made within this framework. 
If assumptions are made, then the conditions under which 
the assumptions are valid can also be made clear. 

Indeed, we have shown in Sec. l7.2l how the general the- 
ory applies to the simple case of non-equilibrium dynamics 
of single-component membranes and under what limiting 
conditions the general theory reduces to the earlier de- 
scriptions |6I7| . Another obvious case to apply the gen- 
eral theory to is the system of active membranes studied 
already both experimentally and theoretically 2 . That 
application will be presented in its full detail in a forth- 
coming paper In the context of this paper, however, 
a question naturally arises concerning that application: 
does the general theory lead to any results that are quali- 
tatively different from the results derived from the theory 
developed in a more intuitive way in Ref. [2]? The answer, 
briefly stated, is yes. The general theory differs conceptu- 
ally from the earlier theory in terms of several new ele- 
ments: first, the inclusions of the bulk viscous stresses in 
the driving force for the transverse permeation processes, 
as stated in Eq. (|1U9|I . as well as in the driving forces 



for lateral diffusion processes, as described in Eq. H138|) : 
secondly, a new proposal for the conditions of matching 
between the bulk and the surface hydrodynamic veloc- 
ities, as formulated in Eq. H11U|) : finally, the prediction 
of a mechanism where surface-tangential gradients of the 
bulk pressures can drive lateral diffusions, as discussed in 
Sec. 17.31 Under the physical conditions that are assumed 
to hold for the active membranes considered, in particu- 
lar, the condition that the membranes may be considered 
impermeable, the first element may be ignored in prac- 
tice. It is much harder to ignore the other two. In fact, 
the analysis in Ref. 35 seems to suggest that the last ele- 
ment may contain a plausible theoretical interpretation of 
certain experimental observations on diffusion dynamics 
of the active proteins in the membranes. 

The general theory may prove to have another poten- 
tial, which is also demonstrated in its application to the 
active membranes 35 . The potential derives from the un- 
derlying philosophy of the theory: that physical effects as- 
sociated with the microscopic "bulk" regions in contact 
with the membrane can be reformulated in terms of sur- 
face excess quantities, if the microscopic details are not the 
focus of a description. The potential is fully taken advan- 
tage of in the application and renders technically tractable 
a calculation that would otherwise be prohibitively diffi- 
cult, if not completely impossible. 

The current theory, despite being sufficiently general, 
is still limited in its scope of coverage of systems or phe- 
nomena. We would like to think that its philosophy will 
lend itself to future development of even more general the- 
ories for describing more complex non-equilibrium phe- 
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A Transforming density fields 

This short appendix is included to describe a technical 
point which may sometimes be encountered in a specific 
calculation based on the theory. In certain situations, it 
can be convenient to do calculations by using different 
linear combinations of the physical density fields, n^'s, 
rather than the fields themselves. An example is that one 
may want to use, as the "working" density fields, those 
combinations that have a certain parity when the chosen 
direction of the normal director to the dividing surface is 
reversed.* 



* See for example Ref. 0, where the so-called "sum field" 
and the "difference field" are used. 
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A general linear transformation of the n^i's takes the 
foUowing form, 

?i'A = ^nBAf^^ , (144) 

B 

where is a constant matrix with det ^ 0. If 

n' AS are chosen to be the "working" density fields, equa- 
tions governing the dynamics of the n'^-fields become the 
equations to be used in the calculation. In order that 
the physics of the problem remain invariant under the 
transformation, other related fields, constants and cur- 
rents must also be transformed. Specifically, the continuity 
equations governing the dynamics of n' as will in fact have 
the same mathematical forms as those for the n^-fields, 
but, different effective masses, m'^'s, must be associated 
with the n'^-fields, which are given by the transformation 

m'^ ^^{M-^)^Bm^ , (145) 

B 

where (M~^)^^ is the inverse of Af i.e. 

B 

Moreover, the mass currents associated with the n'^-fields 
must also be the transform of the physical currents j% 
3 A r' jld ^^'^ ^A,K by the matrix M^^, while the asso- 
ciated chemical potentials are the results of the transfor- 
mation of /z"* by the matrix [M~^)^g. 

The transformations of quantities such as and ki- 
netic constants fiA,B, where the species indices A and B 
do not take value O, are in general more complicated. If 
the following conditions, 

AfO=0, Af4, = 0, (147) 

are satisfied, however, f2A,B still transform according to a 
relatively simple rule 

Q'ab = ^cdM^aM'^b ■ (148) 

It may be worth mentioning that the convention which 
has been used so far of placing the species-label index 
A as either a superscript or a subscript has been chosen 
such that it conforms with the general convention of index 
contraction in matrix algebra. 
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